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Abstract 

In this proceedings, we present STAR's measurement of v-i difference between posi- 
tively charged pions and negatively charged pions at low transverse momentum for Au + 
Au collisions at ^/snn = 200 GeV. The integrated V2 of ir + 0.15 < px < 0.5 GeV/c, 

decreases (increases) linearly with the increasing charge asymmetry A± . The V2 difference 
between n + and tt~ is proportional to A± and the slope parameters have their order of 
magnitude of about 0.01. All of these observed features are consistent with the predica- 
tion based on Chiral Magnetic Wave. The centrality dependence of the V2 difference is 
different than that was predicted. 



1 Introduction 

In relativistic heavy ion collisions, a significant fraction of incident energy is deposited in the 
collision region. Due to the existence of the large energy density, a phase of extremely hot and 
dense matter consisting of quarks and gluons, the Quark Gluon Plasma (QGP) [TJ El El E] , is 
expected to be created. On the other hand, in non-central collisions, spectators keep moving 
along the incident directions. Since spectators are charged objects, moving spectators will 
generate electric current which induces a magnetic field perpendicular to the reaction plane 
(defined by the line connecting the colliding nuclei and the incident direction). It has been 
estimated that the magnetic field created in heavy ion collisions can be as strong as up to 
~ 10 15 Tesla at the beginning of collisions while it decreases rapidly with time [5]. With QGP 
under a such strong magnetic field, it is pointed out that the chiral magnetic effect (CME) 
[21 El El E] and the chiral separation effect (CSE) [TUl E] are expected to exist. The CME in 
Quantum Chromodynamics (QCD) is that the external magnetic field induces a vector current 
with a finite axial chemical potential fi^, which describes the asymmetry between the densities 
of left- and right-handed quarks. 
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If there exist meta-stable Pand CP odd domains, space-time regions occupied by a classical field 
with a nonzero topological charge, under CME it would lead to the asymmetry in the emission 
of positively and negatively charged particles perpendicular to the reaction plane in non-central 
collisions. Such effect has been intensively studied by several experiments [TJl [131 EH E5] ■ The 
other effect in the discussion, CSE, is the separation of chiral charge along the axis of the 
external magnetic field at finite density of vector charge. The resulting axial current is given 
by 

U = ^ V B. (2) 



3v = ^VaB. (1) 
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The local fluctuation of axial current from Eq. (J2]) will induce a local fluctuation of the axial 
chemical potential, and thus according to Eq. (pQ) a fluctuation of electric current. The resulting 
fluctuation of electric charge density will then in turn, again, induce an axial current according 
to Eq. (j2J). With alternating CME and CSE this cycle continues, and it gives rise to a gapless 
excitation regarded as Chiral Magnetic Wave (CMW). Please note that the CMW does not 
require the local parity violation although one of its ingredients, namely the CME, is responsible 
for making the local parity violation signal visible in relativistic heavy ion collisions. 

The CMW has a consequence that also could be observed in experiments. As pointed out in 
a recent theoretical work [16], a CMW in QGP will form a electrical quadrupole moment which 
can lead to more positive charge near the poles of the created fireball and more negative charge 
near the equator. That means, this configuration will cause a difference in elliptic flow, v 2 
[PT] . between positively and negatively charged particles. Taking into account the quadrupole 
moment, the azimuthal angle distribution of positive and negative particles is given by Eq. ([3]), 

« N± [1 + 2v 2 cos(20) =f A±r cos(20)] 

Here N + (AL) is number of positive (negative) particles, A± = (N + — N-)/(N + + AL) is the 
net charge asymmetry and r = 2q e /p e , while q e represents the electrical quadrupole and p e is 
charge density. Authors in [16] also pointed out this difference in elliptic flow should be see via 
7r + and 7r~ due to their similar absorption cross sections in hadronic matter at finite baryon 
density. The elliptic flow difference between tc + and tt~ is expected to be proportional to the 
net charge asymmetry as described in Eq. (j4j). 

Aw 2 CMW = v 2 (tt-) - v 2 (n + ) « rA± (4) 

There are other models that also predicted different elliptic flow between particles and anti- 
particles. Quark transport model predicted v 2 order of particles and anti-particles according 
to their quark components [18] . This model is base on quark coalesce and assumes transported 
quarks have lager v 2 than produced quarks. AMPT model with hadronic potential also pre- 
dicts different v 2 of particles and anti-particles [19]. However, as currently implemented the 
two models mentioned above did not provide a net charge asymmetry dependency of the v 2 
difference. 

In STAR experiment, the difference in integrated v 2 between particles and anti-particles 
has been observed. In this paper, we aim to study such difference differentially as a function 
of net charge asymmetry, and compare that to features predicted by the CMW. 



2 Data Analysis and Results 

STAR detector complex [20] has large acceptance and full azimuthal coverage, which is ideal 
for studying elliptic flow. In this analysis, we used ~ 238M minimum-bias Au + Au events at 
^/s^ = 200 GeV taken by STAR during year 2010. STAR Time Projection Chamber (TPC) 
[2~T] is used to reconstruct tracks and identify particles. All particles used in this analysis are 
required to be in the middle rapidity range < 1.0. For flow analysis, pions are identified 
by ionization energy loss in TPC and required to have 0.15 < px < 0.5 GeV/c. All charged 
particles within 0.15 < px < 12 GeV/c are used to calculate net charge asymmetry except low 
Pt protons, pt < 0.4 GeV/c. Those protons are dominated by protons produced from knock- 
out/nuclear interactions of pions with inner detector material, and thus rejected. Antiprotons 
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within the same pr range are also excluded from the net charge asymmetry calculation in order 
to balance the cut on protons. 

To estimate the elliptic flow of tt + and tc~, we use Q-cumulants method, which is also called 
direct cumulants method [22] ■ Q-cumulants method estimates different order of harmonics via 
multi-particle correlations. This relatively new method needs only on pass over tracks and 
has comprehensive detector inefficiency corrections. In practice of this analysis, two-particle 
correlation is used to estimate t>2{2} of pions and an r/-gap of 0.3 pseudo-rapidity unit on each 
side is applied to suppress the short-range correlation. 

With in a certain centrality bin, we calculate the net charge asymmetry A± event by 
event and divide the sample into five sub-groups so that each sub-group has roughly the same 
number of events. Considering finite detector acceptance and tracking inefficiency, A± obtained 
directly from final state hadrons is different than the true A±, thus it needs to be corrected. 
In discussion below, we call the A± before the correction the observed A±, and otherwise, the 
true A±. The correction procedure on the observed A± will be discussed later. Fig. [1] shows 
the observed A± distribution of 30-40% most central Au + Au events at v /J/vjv = 200 GeV, in 
which red lines divide the statistics equally into five bins. 
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Figure 1: Observed net charge asymmetry A± distribution of 30-40% most central Au + Au 
events at a /J/vjv = 200 GeV. The vertical dashed red lines show the sub-groups of events. 



We measured the integrated t>2{2} of 7r + and n with 0.15 < pt < 0.5 GeV/c and present 



it as a function of observed A±. In Fig. 2(a) , t> 2 of ir + and 7r~ both show a linear relationship 
with respect to the observed A±. With the increasing observed A±, v% of ir~ increases and V2 of 
7r + decreases. As a consequence, the w 2 difference between 7r + and 7r~, Av 2 = V2(tt^) — V2(n + ), 
increases linearly with the increase of the observed A±. All of the features mentioned here are 
consistent with the prediction made in [T6] . 

The slope of the linear relationship between At> 2 and A± is of special interest because 
it has a connection to the electrical quadrupole induced by a CMW in QGP. However, one 
needs to correct the observed A± before a meaningful slope can be extracted. Applying the 
correction on observed A± is conducted in two steps. First, we measure the tracking efficiency 
of charged pions in embedding data and use it as the tracking efficiency of all charged hadrons. 
This efficiency reflects the combined effect of detector and track reconstruction inefficiency. 
Secondly, we generated some Monte-Carlo events by using the event generator Hijing [23]. 
Net charge asymmetry in Hijing events is measured before and after applying the tracking 
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Figure 2: (a) Integrated t>2{2} of n + and n~ , 0.15 < px < 0.5 GeV/c, shows as a function of 
observed A±, 30-40% most central Au + Au events at ^snn = 200 GeV/c. (b) v% difference 
with respect to corrected net charge asymmetry A± fits to a straight line. 



efficiency and all of track selection cuts as applied on real data. The relationship between the 
two measurements of net charge asymmetry is used to map the observed A± to the true A±. 
Roughly, the corrected charge asymmetry A± is a half of its corresponding observed A±. 

In Fig. 2(b) , the At> 2 is plotted a function of A±. A clearly linear relationship is observed. 
Fitting the data points to a straight line, the slope parameter can be extracted. This procedure 
is repeated for all centrality classes. 




% Most Central b (fm) 

Figure 3: (a) Slope extracted from the relationship of At>2 with respect to A± for all centralities. 
Statistical errors only, (b) Slope parameters compare to the theoretical prediction made in [16] , 
which is shown by the solid lines. 



The slope parameters measured for events in the nine centrality bins are shown in Fig. 
Here only statistical uncertainties are shown. It can be seen that slopes have the order of 



magnitude of about 0.01, which is consistent with the theoretical prediction based on CMW. 
The slopes are relatively large in mid-centrality and relatively small in peripheral and central 
collisions. In order to quantitatively compare to the prediction made in [16], we map out 
our centrality definition to impact parameter b according to a previous STAR study based 
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on Monte-Carlo simulation [23]. The comparison is shown in Fig. 3(b) It can be seen 
that, although the magnitude of slopes in general agree with the prediction, the centrality 
dependency is a little different - the data is peaked at a larger b than the prediction. As the 
strong correlation exist between b and the magnetic field, this observation calls for a more 
advanced calculation on the magnetic field that is used to calculate the At> 2 . 

3 Summary 

We measured the v 2 of tt + and ir~ of Au + Au collisions at y / J/v r /v = 200 GeV in differ- 
ent centrality and different net charge asymmetry A± bins. The integrated v 2 of 7r + 
0.15 < pt < 0.5 GeV/c, decreases (increases) linearly with the increasing A±. The v 2 differ- 
ence between ir + and n~ is proportional to A± and the slope parameters have their order of 
magnitude of about 0.01. All of these observed features are consistent with the predication 
based on CMW in [IB]. However, the predicted trend of the slope with respect to impact 
parameter is different than the measurement. 
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